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SUMMAET 


The  tunnel— Induced  upwash  field  for  a  point  element  of  lift  in  a 
rectangular  tunnel  is  shown  to  consist  of  three  superimposed  fields 
that,  except  for  relative  position,  are  independent  of  the  lateral  or 
longitudinal  location  of  the  lift  element  in  the  tunnel.  One  of  these 
fields  is  also  independent  of  the  type  of  tunnel  (open  or  closed)  and 
of  the  vidth-hei^t  ratio  of  the  tunnel,  eind  the  other  two,  which 
depend  on  this  ratio,  are  identical  (except,  perhaps,  for  sign) .  A 
contour  chart  of  the  first  field  is  given  in  the  present  paper;  hence, 
only  one  other  contour  chart  need  he  calculated  for  any  given  tunnel 
to  permit  the  determination  of  the  induced  upwash  field  for  any  position 
of  the  lifting  element.  Contour  charts  of  this  other  field  are  given 
for  three  specific  tunnels:  an  open  tunnel  of  2:1  wldth-hiel^t  ratio, 
a  closed  tunnel  of  2:5  vldth-height  ratio,  and  a  closed  tunnel  of  10:7 
width-hel^t  ratio.  By  superposition  of  results  for  various  locations 
of  the  lifting  element,  the  total  field  may  he  found  for  a  wing  of  any 
plan  form  and  with  any  distrlhutlon  of  lift. 

For  tuimels  that  are  not  rectangular  or  that  cannot  he  considered 
approximately  rectangular,  the  corresponding  procedure  requires  the 
preparation  of  a  chart  for  each  of  several  spanwise  locations  of  the 
lifting  element.  Even  this  procedwe  appears  sin^jler  and  more  generally 
applicable  than  the  calculation  of  induced  upwash  for  a  series  of  wing 
spans  and  sweep  angles. 


UmODUCTION 


As  is  well  known,  the  calculation  of  subsonic  wind-tunnel-wall 
interference  at  a  stral^t  unyawed  lifting  line  is  reducible  to  a 
relatively  simple  two-dimensional  flow  problem;  whereas  the  corresponding 
calcrilatlon  for  a  yawed  or  swept  lifting  line  or  the  calculation  of 
induced  camber  or  of  the  downwash  correction  at  the  tall  cannot  be 
similarly  simplified.  Because  of  the  present  intensive  study  of  swept— 
wing  and  ti’iangular-wing  configurations,  much  effort  is  being  directed 
toward  evaluation  of  tunnel  interference  for  such  wings.  In  general 
the  calculations  are  very  cumbersome  (see,  for  example,  references  1 
and  2) ,  not  only  because  of  the  three-dimensional  character  of  the  flow 
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protlem  iDut  also  'because  the  preparation  of  a  comprehensive  set  of 
corrections  for  a  particular  wind  tunnel  entails  calculations  for  at 
least  a  two—pararaeter  family  of  wings,  that  is,  for  wings  of  various 
spans  and  various  sweep  angles. 

In  some  recent  studies  made  at  the  Langley  Laboratory  it  became 
evident  that  mating  coii^)utations  of  tunnel— Interference  flows  for  a  two— 
parameter  family  of  wings  was  an  unnecessary  complication.  Eesults  of 
equal  accuracy  can  be  achieved  by  a  somewhat  different  and  much  more 
flexible  method  which,  in  general,  requires  computation  of  only  a  one- 
parameter  family  of  charts  and,  for  a  rectangiilar  tunnel,  requires  the 
computation  of  only  two  charts.  Furthermore,  the  coirq)utations  for 
these  charts,  which  can  be  used  for  any  wing,  are  generally  Bliiq)ler 
than  the  computation  of  the  tunnel  Interference  field  for  a  particular 
wing  by  the  usual  method  (in  which  the  wing  loading  is  represented  by  a 
combination  of  yawed  horseshoe  vortices) . 

The  proposed  procedure  is  possibly  Imown  at  some  aeronautical 
laboratories.  Because  of  its  apparent  absence  from  the  literature, 
however,  and  because  procedures  similar  to  those  of  references  1 
and  2  appear  to  be  in  general  use,  the  present  paper  has  been  prepared 
which  outlines  the  method,  describes  the  computations,  and  gives 
examples  of  the  derived  charts.  This  paper  is  concerned  exclusively 
with  the  calculation  of  tunnel— induced  vertical  velocities  in  the 
horizontal  center  plane  of  the  tunnel,  corresponding  to  a  specified 
load  distribution  on  the  wing,  which  is  also  assumed  to  lie  in  the 
same  plane.  Modification  of  the  procediu'e  would  be  required  for 
application  to  triangular  or  highly  swept  wings  at  very  hl^  angles 
of  attack.  No  effort  is  made  to  discuss  the  corrections  to  the 
measured  force  and  moment  characteristics,  since  the  effects  of  the 
Induced  upwash  velocity  on  these  characteristics  are  discussed  in 
reference  1. 
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SYMBOIS 

strength  of  horseshoe  vortex 

span  of  horseshoe  vortex 

lift  of  wing  segment 

area  of  wing  segment 

mean  lift  coefficient  of  wing  segment 

wing  lift  coefficient 
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a,  b,  c,  and  d 

a’,  b*,  aad  d’ 
P,  and  7 


air  density 

stream  velooity 

tunnel  width 

timnel  hei^t 

longitudinal  coordinate 

lateral  coordinate 

vertical  coordinate 

upwash  velocity 

funnel  cross-eectional  area 

cross-sectional  area  of  tunnel  having  same 
proportions  as  actual  tunnel,  but  in  which 
tunnel  width  is  unii^ 

semi— infinite  doublet  line,  similar  to  doublet  line 
in  tunnel 

semi— Infinite  doublet  line,  reverse  of  doublet  line 
in  tunnel 

semi— infinite  source  line 
semi— infinite  sink  line 

points  on  wing  where  lift  is  assumed  to  be 
concentrated 

nearest  lateral  images  of  points  a,  b,  c,  and  d 

points  on  wing  where  tunnel— induced  upwash  angles 
are  to  be  determined 


CALCULATION  OF  TIMEL  INIEEILKENCE 
Representation  of  Wing  Loading 


For  purposes  of  computing  the  tunnel  interference  by  the  method 
to  be  described,  the  assumed  loading  on  the  wing  is  approximated  by 
a  distribution  of  point  concentrations  of  lift,  about  as  Indicated 
in  figure  1.  Roughly,  this  distribution  is  chosen  by  considering  the 
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wing  area  to  te  made  up  of  several  smaller  areas  or  segments^  estimating 
the  lift  on  each^  and  locating  each  equivalent  point  concentration  of 
lift  at  the  approximate  centroid  of  the  lift  on  each  segment.  As  has 
frequently  been  shown,  the  tunnel  interference  field  is  determined  mainly 
hy  the  total  lift  and  the  total  rolling  moment  and  is  otherwise  relatively 
Independent  of  the  precise  lift  dlstrihution  (references  3  aiid  k); 
accordingly,  representing  the  continuous  loading  hy  several  discrete 
point  concentrations  in  the  indicated  manner  is  normally  satisfactory  for' 
the  calculation  of  the  tunnel  Interference  field.  In  any  event,  where  a 
question  arises  as  to  the  adequacy  of  the  representation  used  (as  in  the 
case  of  large— span  wings),  accuracy  can  he  improved  hy  Increasing  the 
number  of  points  or  even  extending  the  procedure  from  a  summation  to  a 
graphical  or  numerical  integration. 

Associated  with  each  concentration  of  lift  is  a  horseshoe  vortex 
of  infinite  strength  and  zero  span  extending  downstream  from  the  point 
where  the  lift  is  considered  to  he  concentrated.  The  moment  PAs  of 
each  horseshoe  vortex  is  given  hy  the  lift  equation  AL  =  pYPAs.  Since 
the  field  of  such  a  degenerate  horseshoe  vortex  is  easily  shown  to  he 
equivalent  to  that  of  a  line  of  source— sink  doublets  (reference  5)^  it 
will  he  referred  to,  for  convenience,  as  a  doublet  line.  The  problem  to 
he  discussed  in  the  succeeding  sections,  then,  is  the  determination  of 
the  tunnel— interference  flow  resulting  from  the  presence  of  a  group  of 
doublet  lines  similar  to  that  indicated  in  figure  1. 


Eectangular  Tunnels 

Image  system  and  interference  field  for  closed  tunnel.—  Figure  2(a) 
shows  the  image  system  for  one  doublet  line  located  in  the  horizontal 
plane  of  symmetry  of  a  closed  rectangular  tunnel.  The  tunnel  boundary 
is  indicated  hy  heavy  lines  and  the  boundaries  of  the  image  tunnels  are 
indicated  hy  light  lines.  The  image  doublets  are  represented  hy  plus  or 
minus  signs  according  as  they  are  the  same  as,  or  the  reverse  of,  the 
doublet  in  the  tunnel.  This  image  system  satisfies  the  boundary  condition 
that  the  velocity  components  normal  to  the  walls  must  he  zero. 

The  complete  system  of  doublets  in  figure  2(a)  is  seen  to  comprise 
two  superimposed,  doubly  infinite,  2b  by  h  rectangular  arrays  of  doublets. 
One  array.  Indicated  by  circles,  may  be  considered  as  centered  at  the 
original  doublet  (double  circle)  in  the  tunnel;  the  other  array,  indicated 
by  squares,  may  be  considered  as  centered  at  the  nearest  horizontal  image 
doublet  (double  sq-uare)*. 

The  Interference  field  is  thus  made  up  of  two  parts: 

(a)  The  field  of  a  complete  rectangular  array  having  its  center 
at  the  double  square. 
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(b)  The  field  of  a  complete  rectangular  array  having  its  center 

at  the  double  circle  with,  however,  the  field  of  the  center 
doublet  omitted  (since  It  represents  the  lifting  element 
Itself) . 

Basis  of  proposed  calctilatlon  procedure.—  The  dimensions  of  each  array, 
and  hence  Its  associated  flow  field,  are  determined  only  by  the  dimensions 
of  the  tunnel  and  are  Independent  of  the  lateral  or  longitudinal  location 
of  the  lifting  element  In  the  tunnel.  Accordingly,  once  the  two  fields 
have  been  calculated,  they  may  be  used  for  finding  the  Interference 
corresponding  to  a  lifting  element  located  anywhere  In  the  horizontal 
center  plane  of  the  given  closed  rectangular  tunnel.  The  procedure  is 
Indicated  in  figure  3  which  shows  (plan  view)  a  lifting  element  and 
Its, nearest  image.  In  figure  3(a)  is  indicated  the  contour  chart  of 
induced  upwash  velocities  calculated  for  a  2b  by  h  doubly  Infinite 
array  of  unit  doublets  with  the  center  doublet  omitted.  The  point  on 
the  contour  chart  that  is  located  at  the  head  of  the  omitted  doublet  line 
is  indicated  as  the  origin;  and  the  chart  is  placed  so  that  this  point  falls 
on  the  lifting  element,  designated  a.  In  figure  3(b)  is  indicated  the 
contour  chart  for  the  complete  doubly  infinite  array,  placed  so  that  its 
origin  falls  on  the  first  image,  designated  a’.  At  any  specified  point  a 
in  the  horizontal  center  plane  of  the  tunnel,  the  induced  upwash  corre¬ 
sponding  to  the  given  element  of  lift  is  found  by  adding  the  values  read 
at  that  point  from  the  contour  charts  in  figures  3(a)  and  3(*t))  and 
multiplying  the  s^um  by  the  strength  of  the  equivalent  doublet. 

The  procedure  may  be  slightly  modified  to  take  advantage  of  the  fact 
that  the  chart  of  figure  3(b)  is  equivalent  to  the  sum  of  two  other  charts, 
namely,  the  chart  of  figure  3(0*)  ^  which  is  for  the  doubly  infinite  array 
with  the  center  doublet  omitted,  and  a  chart  for  a  single  doublet  (that 
which  is  omitted  in  the  chart  of  figure  3(0')).  Accordingly,  the  chart  of  • 
figure  3(a)  and  a  chart  for  a  single  doublet  should  suffice  to  obtain  the 
desired  upwash  values.  In  this  modification,  the  step  indicated  in 
figure  3(b)  is  replaced  by  the  two  steps  indicated  in  figures  4(a)  and  4(b). 
Three  readings  are  thus  necessary  instead  of  two;  since  the  chart  for  a 
single  doublet  is  given  in  the  present  paper  (fig.  5)y  however,  this 
modification  requires  that  only  the  chart  of  figure  3(0)  be  prepared  for 
each  given  rectangular  tunnel. 

Image  systems  for  other  rectangular  tunnels.—  In  figures  2(b)  to  2(f), 
are  shown  the  image  systems  for  five 'other  rectangular  tunnels,  namely, 
those  that  are,  respectively, 

(1)  Open  on  all  four  sides 

(2)  Closed  at  the  sides  but  open  at  the  top  and  bottom 

(3)  Closed  at  the  top  and  bottom  but  open  at  the  sides 

(4)  Closed  only  at  the  bottom 

(5)  Closed,  containing  a  semlspan  reflection  model 
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Open  "boundaries  are  indicated  "by  dashed  lines  and  closed  houndarles,  "by 
solid  lines.  In  each  of  these  cases,  as  with  the  completely  closed 
tunnel,  the  image  aysteni  is  seen  to  he  composed  of  two  2h  hy  h  rectangular 
arrays  so  that  the  procedure  Just  outlined  should  also  apply  for  these 
configurations,  except  that  for  some  of  the  tunnels,  however,  the  two 
arrays  have  opposite  signs  so  that  the  readings  from  figures  4(a) 
and  4(h)  must  he  subtracted  from,  instead  of  added  to,  those  from 
figure  3(a) .  The  rectangular  arrays  are  not,  however,  all  of  the  same 
type  as  for  the  conpletely  closed  tunnel.  The  image  systems  can  he 
briefly  described  as  follows: 

For  a  closed  timnel  (fig.  2(a)),  each  array  consists  of  alternate 
horizontal  rows  of  plus  and  minus  doublets,  and  the  two  arrays  have  the 
same  sign  (that  is,  the  nearest  horizontal  image  is  plus). 

For  an  open  tunnel  (fig.  2(h)),  all  the  doublets  of  each  array  have 
the  same  sign,  hut  the  two  arrays  have  opposite  signs. 

For  a  tunnel  closed  at  the  sides  hut  open  at  the  top  and  bottom 
(fig.  2(c)),  all  the  doublets  of  each  array  have  the  same  sign  and  the 
two  arrays  also  have  the  same  sign;  that  is,  only  plus  doublets  occur 
in  the  image  system. 

For  a  tunnel  closed  at  the  top  and  bottom  but  open  at  the  sides 
(fig.  2(d)),  each  array  consists  of  alternate  rows  of  plus  and  minus 
doublets,  and  the  two  arrays  have  opposite  signs. 

For  a  tunnel  closed  only  at  the  bottom  (fig.  2(e)),  the  signs  for 
the  image  system  in  the  first,  third,  fifth,  .  .  .  rows  above  the  plane 
of  the  wing  are,  respectively,  opposite  to  those  in  the  first,  third, 
fifth,  .  .  .  rows  below  the  plane  of  the  wing.  Since  the  net  effect 
of  the  odd  numbered  rows  on  the  upwash  velocity  in  the  center  plane  is 
thus  zero,  they  may  be  neglected.  Two  2b  by  2h  arrays  remain  having 
alternate  rows  of  plus  and  minus  doublets  and  the  two  arrays  have 
opposite  signs. 

For  a  closed  txmnel  containing  a  semispan  reflection  model 
(fig.  2(f)),  each  array  has  alternate  rows  of  plus  and  minus  doublets 
and  the  two  arrays  have  the  same  sign.  In  the  indicated  image  system 
the  image  in  the  double  square,  which  represents  the  second  half  of 
the  reflection  model,  must,  like  the  model  Itself,  be  excluded  from 
consideration  in  calculating  the  interference  field.  The  reflection 
system  thus  consists  of  two  arrays  in  each  of  which  the  central  doublet 
is  omitted.  For  this  case,  then,  only  one  chart  (similar  to  that  of 
fig.  3(a))  is  used.  The  same  conclusion  applies  whenever  a  semispan 
reflection  model  is  tested  with  a  tunnel  wall  as  the  reflection  plane. 

As  a  corollary,  it  follows  that,  if  a  tunnel  is  used  to  test  only 
symmetrical,  unyawed,  full— span  models,  only  one  contour  chart  need  be 
prepared,  and  the  chart  would  represent  the  field  of  a  b  by  h  (instead 
of  a  2b  by  h)  array.  This  conclusion  can  be  readily  verified  by 
examining  figure  2(g)  which  represents  a  closed  timnel  containing  a 
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full— span  aymmetrical  model  and  viiich  is  Identical  with  figure  2(f)  except 
that  the  reflection  plane  and  all  its  images  are  omitted. 

Basic  formulas  and  summation  procedures  for  calculating  contour  charts. 
The  potential  of  a  unit  doublet  at  the  origin  with  its  axis  vertical  is 


4it(x2  +  72  +  5.2)3/2 

The  potential  of  a  doublet  line  extending  along  the  x-axls  from  the  origin 
to  infinity  is  then 
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The  corresponding  upwash  velocity  is 
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For  any  particular  array  of  unit  doullets^  the  upwash  velocity  at  a 
point  (the  contour  value  at  that  point  on  the  contour  chart)  is  the  sum 
of  the  values  given  hy  this  formula  for  a  series  of  values  of  y  differing 
hy  213  and  a  series  of  values  of  z  differing  hy  h,  with  appropriate 
signs  according  to  the  type  of  array.  Such  a  double  series  usually  con¬ 
verges  rather  slowly,  however,  and,  in  general,  the  practicability  of  the 
summation  depends  on  the  use  of  certain  approximation  methods  for  summing 
the  fields  of  all  but  an  inner  group  of  doublets  surrounding  the  origin. 
These  approximation  methods,  which  are  very  similar  to  those  used  in  the 
two-dimensional  studies  of  reference  6,  are  reviewed  in  the  following 
paragraphs . 

The  field  of  a  doublet  line  is  approximately  the  same  as  the  field 
of  a  horseshoe  vortex  of  the  same  moment,  provided  the  distance  from  the 
doublet  line  to  the  point  where  the  field  is  being  considered  is  suf¬ 
ficiently  large  relative  to  the  span  of  the  horseshoe  vortex.  Thus,  in 
computing  the  field  near  y  =  0,  a  row  of  equal  doublets  at,  for  Instance, 
y  =  4b,  6b,  8b,  .  .  .  may  be  replaced  by  a  row  of  horseshoe  vortices  of 
span  2b  having  their  trailing  vortices  at  y  =  3b  and  5^^  51^  and  Jh, 

7b  and  9b,  ...  In  this  representation,  all  the  trailing  vortices  except 
the  innermost  one  cancel  in  pairs,  so  that  the  infinite  row  of  doublets  is 
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eq.ulTalent  to  an  L-vortex  of  vhlcli  the  trailing  portion  is  at  y  =  31> 
and  of  which  the  hound  portion  extends  along  the  y-axis  from  3h  to 
infinity.  The  field  of  this  infinite  Ir^ortex  is  easily  calculated  hy 
the  Biot-Savart  law. 

For  a  horizontal  row  of  doublets  lying  a  reasonable  distance  above 
the  origin,  all  may  be  replaced  by  horseshoe  vortices  of  span  2b,  in 
which  case  all  the  trailing  vortices  cancel  in  pairs  and  only  the  bound 
vortex  extending  from  -»  to  +»  remains.  Its  field  is  merely  that  of 
a  two-dimensional  vortex. 

By  means  of  the  approximate  representations  Just  described,  a 
rectangular  array  of  doublets  in  which  the  alternate  rows  have  plus  and 
minus  signs  can  be  assumed  approximately  equivalent  to  an  inner  group  of 
doublets  around  the  origin  (those  that  are  too  close  to  be  adequately 
replaced  by  horseshoe  vortices)  and  an  outer  arrangement  of  L-vortices 
and  two-dimensional  bound  vortices  (fig.  6(a)).  Because  of  the  alter¬ 
nating  signs,  the  upwash  velocity  at  any  point  due  to  the  two-dimensional 
boiuid  vortices  may  be  formulated'  as  the  sum  of  the  terms  of  an  alternating 
series,  which  can  be  readily  eval'uated. 

For  a  rectangular  array  of  doublets  in  which  all  have  the  same  sign, 
a  different  replacement  system  is  more  convenient.  Instead  of  being 
extended  horizontally  into  a  horseshoe  vortex,  the  doublet  is  extended 
vertically  into  a  source  line  and  a  sink:  line,  a  distance  h  apart. 

The  source  lines  and  sink  lines  in  any  column  cancel  each  other  in 

pairs  and  only  the  source  or  sink  lines,  at  a  distance  -g  above  or  below 

the  inner  group  of  doublets,  remain.  This  equivalent  representation  for 
a  rectangular  array  in  which  all  doublets  have  the  same  sign  is  shown  in 
figure  6(b) . 

It  may  be  noted  that  the  previously  mentioned  L-vortices  are  acfually 
horseshoe  vortices  of  which  the  second  trailing  vortex  Is  at  y  =  •  and 
that  similar  trailing  vortices  at  y  =  are  associated  with  the  infinite 
bound  vortices.  Similarly,  where  the  source  lines  and  sink  lines  in  the 
columns  cancel  each  other  in  pairs,  a  row  of  sink  lines  remain  at  z  -  +» 
and  a  row  of  source  lines  at  z  =  — »  The  trailing  vortices  at  y  =  •±-“ 
require  no  special  consideration  because  their  alternating  signs  would, 
in  any  event,  resiilt  in  zero  net  effect  in  the  region  of  the  origin.  The 
sink  and  source  lines,  however,  would  be  expected  to  result  in  a  uniform 
upwash  throughout  the  field,  provided  the  extent  of  the  rows  of  source 
and  sink  lines  is  of  a  higher  order  of  inf ini^  than  is  their  dlst^ce 
from  the  origin.  For  the  open  tunnel  (fig.  2(b)),  whether  this  uniform 
vertical  velocity  is  included  or  not  is  Immaterial,  since  the  two  arrays 
have  opposite  signs  and  the  uniform  upwash  would  thereby  be  eliminated  in 
any  case.  The  question  still  remains,  however,  for  the  tunnel  of 
figure  2(c),  where  the  two  arrays  have  the  same  sign.  Actually,  a 
discussion  of  the  relative  orders  of  the  two  infinities  is  not  necessary 
or  desirable  inasmuch  as  a  simple  physical  criterion  is  available  for 
such  cases,  namely,  that  the  total  Induced  upwash  should _ approacn  zero 
at  a  large  distance  upstream  from  the  wing.  Since  the  field  indicated 
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In  figure  6(1))  already  Batiafiea  this  upstream  boundary  condition,  the 
imlform  upwaah  that  wuld  result  from  the  source  and  slrOc  lines  at 
infinity  should  clearly  not  he  Included. 

Eemarks  on  preparation  of  charts  and  scale  factors.—  In  the 
calculations  for  the  charts^  it  is  convenient  to  consider  the  tunnel 
vidth  as  unity  and  to  assume  unit  doublets  in  the  array  (that  is^  to 
use  equation  (l)  for  the  field  of  each  doublet) .  The  half-vidth  of 
the  charts  need  never  exceed  1.5  tunnel  widths^  as  is  obvious  from 
figure  h.  An  example  of  a  complete  charts  used  in  estimating  corrections 
for  the  Langley  full-scale  tunnel  (idealized  to  a  60—  by  3O— foot  rec¬ 
tangular  tunnel),  is  shown  in  figure  7(a)*  The  tunnel  was  considered  as 
a  1  unit  by  l/2  unit  rectangle,  so  that  the  tunnel  area  D  was  l/2 
and  the  doublet  Images  formed  a  2  by  l/2  array.  The  chart  half-width 
is  considerably  less  than, 1.5  tunnel  widths  because  results  for  this 
case  showed  that  the  outer  contour  values,  when  added  to  those  of 
figure  5^  were  insignificant.  In  figure  7(b)  is  shown  a  chart  for  a 
closed  3“”  by  7*5—foot  rectangular  tunnel.  The  chart  was  computed  for 
a  2  by  2.5  array  of  unit  doublets  corresponding  to  a  tunnel  area  D 
of  2.5.  Figure  7(c)  shows  a  chart  similarly  developed  for  7—  by  10— foot 
rectangular  tunnels . 

The  chart  of  figure  5  for  the  single  unit  doublet,  as  computed 
from  equation  (l)  with  z  =  0,  was  plotted  with  the  unit  distance 
indicated  on  the  coordinate  axes.  In  keeping  with  the  preceding 
discussion,  the  chart  half-width  is  1.5  units. 

The  scale  factor  for  the  eventual  application  of  the  charts  is 
developed  as  follows: 

The  lift  AL  associated  with  a  horseshoe  vortex  of  strength  F 
and  sp^  As  (actual  wing  dimensions)  is  given  by 


AL  =  ctAS^I^ 
=  (rAs)pV 


from  which 


VA^.SJfL 


which  is  the  equivalent  doublet  strength  of  a  segment  of  wing  area  AS 
having  a  mean  lift  coefficient  07 .  The  chart  contours  give  the  upwash 
velocity  w  for  a  doublet  of  unit  strength  in  a  tunnel  of  cross- 
sectional  area  D.  For  the  case  of  a  doublet  of  strength  FAs  in  a 
tunnel  of  area  C, 


w  =  Chart  reading  X  FAs  X  ^ 

0 
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or,  vitli  PAs  replaced  lay  c^AS? /2,  the  upwash  angle  in  radians  is 

*  =  Chart  reading  x 


or  in  degrees. 


28. 6c  2^ 
Chart  reading  x  - s - 


If  the  suggested  conrention  —  that  the  tunnel  width  he  taken  as  unity 
in  computing  the  charts  —  is  followed,  then  the  factor  D/C  la 
merely  l/b^  and  the  expression  for  the  upwash  angle  in  degrees  is 

.  28.6c7AS 

Chart  reading  x - 1 — 

h^ 


Use  of  charts  in  computation  of  tunnel  interference  flow.—  Assume 
that  corrections  are  desired  for  a  closed  rectangular  tunnel,  for  which 
a  chart  of  the  type  shown  in  figure  7  (that  is,  for  an  Infinite  array  of 
doublets  with  the  center  doublet  omitted)  has  been  prepared.  This  chart 
is  designated  chart  A.  The  chart  of  figure  5,  for  a  single  doublet,  is 
designated  chart  B.  The  procedure  may  then  be  outlined  as  follows: 

(1)  Sketch  the  complete  wing  to  the  scale  of  the  charts.  Show  the 
sides  of  the  tunnel  and  the  nearest  Images  of  the  semispans  (fig.  8). 

(2)  Assume  the  lift  to  be  concentrated  at,  say,  two  points  on  each 
semispan,  and  estimate  the  lift  at  each  point  in  terms  of  c^AS.  The 
sum  of  the  foiu*  values  of  C2AS  must  eq.\:ial  Cxp  for  the  complete  wing. 
Also,  if  rolling  moments  are  being  considered,  the  rolling-moment 
coefficient  of  the  approximate  representation  shoTild  equal  the  rolling- 
moment  coefficient  of  the  wing.  Indicate  on  the  sketch  the  four  points 
b,  c,  and  d  and  their  nearest  Images  a*,  b',  c',  and  d' .  Also, 
locate  on  the  sketch  the  points  at  which  the  tunnel— Induced  angle  is 
desired,  say  a,  3,  and  7.  In  general,  where  only  a  few  points  are 
used,  greatest  accuracy  is  achieved  by  assuming  the  lift  (in  the  case  of 

an  unflapped  wing)  to  be  concentrated  along  the  chord  line  and 
determining  the  upwash  angles  along  the  chord  line,  as  indicated 
in  figure  8.  (The  drag  correction,  however,  is  probably  more  correctly 
determined  from  the  upwash  angles  at  the  i  — chord  line  where  the  lift  is 
concentrated . ) 

(3)  Place  the  origin  of  chart  A  at  point  a  and  read  the  chart 
contours  at  a,  P,  and  7.  Bepeat  for  points  b,  c,  and  d.  Then  place 
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the  origin  of  chart  A  at  point  a*  and  read  the  chart  contours  at 
a,  P,  and  7.  Eepeat  for  points  h',  c’,  and  d’ .  finally,  place  the 
origin  of  chart  B  at  point  a’  and  read  the  cheart  at  ct,  p,  and  7. 
Bepeat  for  points  h’,  c',  and  d* . 


(h)  Finally,  for  the  upwash  angle  at,  say,  point  a,  find  the 
algehralc  sum  of  the  three  chairt  readings  (for  eiaa^jle,  (a, a) 
and  (a', a)  from  chart  A  and  (aSot)  from  chart  B)  for  each  of  the 
four  points  a,  h,  c,  and  d.  Multiply  each  of  these  four  sums  hy  the 

value  of  — -j;^C2AS  for  the  wing  segment  under  consideration 

(a^  Cy  or  d)  and  add.  The  total  is  the  funnel— induced  upwash  angle^ 
in  degrees^  at  point  a.  Proceed  similarly  for  points  P  and  7. 


When  the  vertical  eymmetry  plane  of  the  tunnel  is  also  the  vertical 
symmetry  plane  of  the  wing  (in  general^  for  zero  yaw  and  zero  aileron 
deflection),  the  work  can  he  somewhat  reduced.  Only  two  load  points  a 
and  h  on  the  wing  and  their  two  reflections  a’  and  h*  need  he  used, 
althou^  the  chart  readings  must  -  still  he  obtained  at  the  two  symmetri¬ 
cally  located  points  a  and  7  on  hoth  sides  of  the  wing.  The  net 
induced  angle  at  a  or  7  is  obtained  hy  adding  the  results  for  a 
and  75  the  net  induced  angle  at  P  is  obtained  hy  doubling  the  result 
f  or  P . 


Nonrectangular  Tunnels 

For  a  circular  or  other  nonrectangular  tunnel,  the  problem  cannot 
he  reduced  to  that  of  calcxilatlng  a  single  chart.  A  series  of  contour 
charts,  giving  the  tunnel— induced  upwash  in  the  horizontal  center  plane 
of  the  tunnel,  must  he  constructed  for  a  series  of  spanwlse  locations 
of  the  doublet  in  the  txinnel.  Even  so,  only  a  single  parameter  —  the 
spanwlse  location  of  the  doublet  —  is  involved,  since  variations  in 
longitudinal  location  of  the  lifting  element  (such  as  those  for  the 
different  lifting  elements  along  a  swept  wing)  are  readily  taken  into 
account  by  longitudinal  shift  of  the  contour  chart.  Erom  a  study  of 
reference  2  it  appears  reasonably  certain  not  only  that  the  calculations 
herein  proposed  would  have  been  easier  than  the  calculations  therein 
Indicated  but  also  that  the  eventual  application  of  the  results  to  wings 
of  irregular  load  distributions  would  also  have  been  easier.  The 
procedure  Indicated  in  reference  7  seems  to  be  of  this  type. 

With  regard  to  studies  similar  to  those  of  reference  2,  it  may  be 
of  value  to  point  out  that,  at  least  theoretically,  a  aeries  of  calcu¬ 
lations  for  any  one  sweep  angle  should  suffice  for  computing  the  tunnel 
effect  for  any  other  sweep  angle.  Figure  9(a)  shows,  for  example,  how 
the  loading  on  a  60^  swept  wing  may  be  approximated  by  a  single  horseshoe 
vortex  and  two  pairs  of  unswept  horseshoe  vortices,  where  the  inner 
vortex  of  each  pair  has  the  same  strength  as  the  superimposed  outer 
vortex  but  has  opposite  rotation.  Figure  9(’t))  shows,  similarly,  how 
a  pair  of  horseshoe  vortices  and  a  single  horseshoe  vortex,  all 
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swept  45°,  might  he  used  to  approximate  the  loading  on  a  60°  swept  wing 
with  sufficient  accuracy  to  calculate  the  timnel  interference  velocities. 


CONCLUDING  REMARKS 


For  a  lifting  wing  of  arhitratry  loading  and  plan  form,  situated  in 
the  horizontal  center  plane  of  a  recteuigular  tunnel,  the  tunnel— Induced 
upwash  velocities  in  the  same  plane  can  he  readily  calculated  with  the 
aid  of  two  charts.  One  is  given  in  th®  present  paper;  the  other  must  he 
computed  for  each  tunnel.  Such  simplification  is  not  possible  for 
nonrectangular  tunnels.  In  any  case,  however,  computations  for  a  series 
of  wings  of  various  spans  and  various  sweep  angles  are  unnecessary.  Thus, 
computations  for  several  spanwlse  positions  of  the  lifting  element  or 
computations  for  a  set  of  unswept  wings  of  various  spans  sho^lld  provide 
a  haais  for  computing  tunnel— Induced  upwash  velocities  for  any  wing  in 
the  horizontal  center  plane  of  the  tunnel. 


Langley  Aeronautical  Laboratory 

National  Advisory  Committee  for  Aeronautics 
Langley  Field,  Va.,  July  27,  19W 
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Figure  1.-  Field  of  a  lifting  wing  represented  as  the  field  of  four  doublet  lines 
extending  downstream  from  four  points  where  the  lift  is  assumed  to  be 
concentrated. 


Figure  2.-  Image  configurations  for  doublet  lines  in  rectangular  tunnels. 
Doublet  in  tunnel  is  indicated  by  double  circle;  nearest  lateral  image  is 
indicated  by  double  square. 
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(b)  Open  tunnel. 
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Figure  2.-  Continued 
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Figure  2.-  Continued. 
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Closed  tunnel  containing  a  full-span 
symmetrical  model. 


Figure  2.-  Concluded. 
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(a)  Chart  representing  field  of  doubly  infinite 

array  of  unit  doublets  with  center  doublet 
omitted;  origin  of  chart  is  placed  on  point  a. 


(b)  Chart  representing  field  of  complete  doubly  infinite 
array  of  unit  doublets;  origin  of  chart  is  placed 
on  nearest  image  a' . 


Figure  3.  -  Timnel  interference  velocity  at  point  a ,  corresponding  to  lift 
concentration  at  point  a,  determined  by  two  contour  charts. 


(b)  Chart  representing  field  of  a  single 
unit  doublet. 

Figure  4.-  Two  contour  charts,  the  sum  of  which  equals  the  chart  of 

figure  3(b). 
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;y  due  to  a  unit  positive  semi -infinite  doublet  line, 
ad  of  the  doublet  line. 


(a)  Array  having  alternate  rows  of  plus  and 
minus  doublets. 


(b)  Array  having  only  plus  doublets. 


Figure  6.-  Approximate  representations  of  the  two  types  of  doubly  infinite, 
rectangular  arrays  of  doublets  with  center  doublet  omitted. 


Plow  direction 


doubly  infinite  array  of  unit  semi -infinite  doublet  lines  with  the  center  doublet  line  omitted.  The 
origin  0  is  the  head  of  the  center  doublet  line. 


Flow  direction 


Langley  two-dimensional  tunnels  (3  by  7.5  ft) ;  calculated  for  a  2  (horizontal)  oy  2.5  (vertical) 
The  doublets  alternate  in  sign  vertically,  with  those  in  the  horizontal  center  plane  positive. 


Flow  direction 
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Figure  8.-  Sketch  used  with  charts  A  and  B  for  determining  tunnel -induced 
upwash  angles.  Sketch  shows  points  a,  b,  c,  and  d  where  lift  is  assiuned 
to  be  concentrated,  the  nearest  images  a',  b',  c*,  and  d',  and  points  a, 
p ,  and  y  where  the  tunnel -induced  upwash  angles  are  to  be  determined. 
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Figure  9. 


(b)  Approximation  by  horseshoe  vortices  of 
45°  sweep. 

-  Approximation  of  the  loading  on  a  60°  sweptback  wing  by  means 
of  horseshoe  vertices  of  different  sweep. 


